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Abstract 

A natural geometry, arising from the embedding into a Hilbert space of the 
parametrised probabihty measure for a given lattice model, is used to study 
the symmetry properties of real-space renormalisation group (RG) flow. In 
the projective state space this flow is shown to have two contributions: a 
gradient term, which generates a projective automorphism of the state space 
for each given length scale; and an explicit correction. We then argue that 
this structure implies the absence of any symmetry of a geodesic type for 
the RG flow when restricted to the parameter space submanifold of the state 
space. This is demonstrated explicitly via a study of the one dimensional 
Ising model in an external field. In this example we construct exact expres- 
sions for the beta functions associated with the flow induced by infinitesimal 
rescaling. These constitute a generating vector field for RG diffeomorphisms 
on the parameter space manifold, and we analyse the symmetry properties of 
this transformation. The results indicate an approximate conformal Killing 
symmetry near the critical point, but no generic symmetry of the RG flow 
globally on the parameter space. 
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I. INTRODUCTION 



The construction of effective tfieories for the low energy modes from their underlying 
microscopic dynamics or, equivalently, the study of the thermodynamic behaviour of spin and 
lattice field systems is a fundamental problem in modern physics. The standard theoretical 
framework for this procedure, the renormalisation group provides a methodology for 
systematically integrating out high momentum or short distance modes leaving an effective 
dynamical system for the long wavelength physics. This effective theory then determines, 
in particular, the critical behaviour of field theories or alternatively the thermodynamic 
behaviour of lattice or statistical mechanical spin systems. 

The fundamental obstacle in following this procedure through to completion is that 
coarse graining the system, regardless of the detailed approach employed, generically leads 
to a highly complex effective theory introducing an often arbitrarily high number of addi- 
tional interactions. While some simplifications occur near the upper critical dimension for 
a given system, allowing systematic techniques such as the e-expansion to be employed, 
generically one requires drastic approximations in order to render the RG procedure prac- 
tically implement able. In particular, we focus on real-space renormalisation [Q] methods 
which implement an intuitive blocking picture in attempting to explain certain universal 
aspects of critical behaviour, often near the lower critical dimension. However, except in a 
few tractable cases, such calculational schemes involve significant truncations, and there is 
rarely a useful criterion for selecting a particular approximation. This is also a significant 
problem in more recent attempts to study field theories via Wilsonian RG techniques as 
the dominant infrared degrees of freedom in such cases may differ greatly from the relevant 
microscopic fields. 

With these problems in mind it is clearly of interest to study the structure of the RG in 
some detail. Generically the Hamiltonian for a spin system, or the Wilsonian action for a field 
theory, at a given scale may be represented in the form, H = Y^iO^Hi. This is generally an 
infinite series of interaction Hamiltonians Hi and their associated scale dependent couplings 
9^{t). The RG acts on these couplings as one coarse grains the system, and one may gain 
insight into this process by studying the geometry of the parameter space, whose singularity 
structure and possible symmetries determine the properties of the RG flow, in particular 
the fixed points. However, the manifold M. of the parameter set dealt with here is typically 
a subspace of an infinite dimensional Hilbert space Ti 0] and this infinite dimensionality, 
which in effect is the cause of the singularities, gives rise to various technical problems that 
have to be treated with caution. Nevertheless, the ability to formulate certain RG related 
questions in a geometric formalism allows the utilisation of various powerful techniques in 
differential geometry. 

The Riemannian structure of this parameter manifold M., as we discuss in more detail 
below, is inherited from the structure of the probability measure of the system over its 
configuration space 0]. Various authors have utilised this structure to study aspects of 
statistical mechanical systems and field theories, and it has been shown that |0 the RG acts 
via diffeomorphisms of this parameter space manifold generated by the /3-function vector 
field. It is then a clear corollary that the geometry can place restrictions on the /3-function 
of the theory if other symmetries of the parameter space exist, or indeed if the RG flow itself 
is a symmetry of M. . 
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In this paper, we investigate this last question, and consider whether the RG diffeo- 
morphism actually corresponds to a symmetry of the parameter space. This is motivated 
by earlier work |^,^ which suggests that an approximate symmetry holds near the critical 
points, and also to some extent by recent work on the geometric structure of the 



state space of statistical mechanical systems, and in particular the parametric evolution of 
the states. A consequence of this latter work is that for a system whose action is express- 
ible in the form Y^iO^Hi, the parametric evolution of the state corresponds to a projective 
automorphism of the projective Hilbert space, which is the most general transformation of 
a Riemannian manifold mapping geodesies into geodesies. 

In generalising this result to systems where all the parameters depend on a single dimen- 
sionful scale, we show in Section 2 that the generating vector field of the flow, when restricted 
to a constant scale, still generates a projective automorphism of the state-space. However, 
in considering the full scale dependence of the states, it is shown that the integral curve 
corresponding to RG evolution only generates a projective automorphism of the state-space 
up to an anomalous correction. Since the parameter space can be viewed as a submanifold 
of the projective state-space, this suggests the absence of such a geodesic type symmetry of 
the RG flow on the parameter space itself. The generality of this construction allows us to 
conjecture that in general the RG does not generate any standard global symmetry of this 
type on the parameter space. 

In order to test this conjecture we provide an explicit counter-example in which both the 
RG trajectory, and the underlying parameter space geometry, can be exactly determined. 
While the RG procedure generally requires various approximations, there are a number of 
exactly tractable models known in statistical mechanics [0,^, and in the present paper we 
study one such model exhaustively by use of real-space RG methods, in order to determine 
the properties of RG flow in the relevant parameter space. 

The example we consider here is the one dimensional Ising model in an external mag- 
netic field. As we shall discuss in more detail below, the parameter space in this case is a 
two dimensional manifold endowed with a Riemannian structure in terms of the Fisher-Rao 
metric. In Section 3 we consider the exact parameter space geometry of this system both 
for a finite A^-spin chain and in the thermodynamic limit. In Section 4 we utilise a transfer 
matrix method to obtain an exact expression for the vector field associated with an infinites- 
imal change of the lattice spacing, the components being the differential RG /^-functions. 
Using this result we study, in Section 5, the symmetry associated with the flow. We find, in 
particular, that the generating vector field corresponds to a conformal Killing vector field in 
the vicinity of the critical point. However, this vector field does not generate such a trans- 
formation, nor the most general additional class of mappings, a projective automorphism, 
globally on the manifold. Therefore we are led to conclude, via this explicit counter-example, 
that there is no generic global symmetry structure of a standard form for the RG flow. We 
finish in Section 6 with some concluding remarks. 



II. STATE SPACE RENORMALISATION 

We begin by considering the effect of renormalisation group flow in the actual state space 
of a given system, which can be viewed as the space of rays through the origin of a real 
Hilbert space Ti, that is, the real projective ra-space i?P" (possibly infinite dimensional). 
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Before studying the behaviour of RG flow in this space we review briefly the notion of the 



state space in statistical mechanics [|TT 



In a purely statistical mechanical context we are given the parametrised family of prob- 
ability distributions, taking the form of the Gibbs measure, 



p{x, 6) = q{x) exp 



(1) 



where the variable x ranges over the configuration space, Hj{x) represents the form of the 
energy, Wg is a normalisation factor, and q{x) determines the distribution at 6^ = 0. In a 
field theoretic context the operators Hj effectively determine the form of the action, and the 
parameters 6^ are viewed as the coupling constants. Now, by taking the square-root of the 
distribution function (|1]) we can map, for each given value of 6\ the space of probabilities 
onto a manifold in a real Hilbert space Ti. Since the probability distribution is normalised, 
we find that for each fixed value of 6^ , the probability state corresponds to a point on the 
unit sphere S dTi. Then, by choosing a suitable basis of vectors in 7i, we can consider the 



vector ijj"'{6) as representing the point characterised by yp{x,6). In other words, we regard 
iIj^{6) as a state vector corresponding to the distribution 

First we would like to formulate a Hilbert space characterisation of this distribution. In 
fact, it can be shown p!0| , p!l|] that the state vector ip'^{9) in TC corresponding to the Gibbs 
distribution ([^) satisfies the differential equation. 



(2) 



where the operators Hjab in H represent the energy Hj{x), Hjab = Hjab — gabE^[Hj], with 



licic Liic ujjciauuis J-J-jab iJ^J- ' L J-cpicsciiL iiic ciicigj' iijyxj, Hjab — -tijab ^u,u~ifji 

E[-] denoting the expectation. The solution to this equation is given by an exponential 
family |^ of states 



^"(0) = exp 



(3) 



where Wq = Wg — Wq and = ?/'"(0) is the prescribed state at 9^ = 0. 

The real Hilbert space Ti, however, is not what we view here as the true state space of 
the physical model under study since there is an extra degree of freedom, namely, the overall 
normalisation. That is to say, the expectation of any physical observable is independent of 
the value of the normalisation. Hence we can gauge this extra degree of freedom away by 
identifying all the points in Ti along the given ray that passes through the origin of Ti., 
corresponding to different normalisation factors. The resulting space obtained is the real 
projective n-space RP^ which we view as the actual state space, with the state vectors 
ip"- in Ti. also representing the homogeneous coordinates for i?P". Note that we encounter 
an analogous situation in quantum mechanics where the overall phase degree of freedom 



associated with a given quantum state can be eliminated ||T2| to recover the quantum phase 
space CP". 

Now, given the exponential states (^, we may wish to ask if there is any symmetry 
associated with the flow induced by changing the canonical parameters 6^ in RP"'. This 
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line of enquiry has been investigated recently |TT|, and it was shown that the induced flow 



gives rise to a projective automorphism on the state manifold — a general transformation 
of the kind that projects geodesies onto geodesies. Furthermore, it has also been shown 
that the induced flow is a Hamiltonian gradient flow with respect to the natural 'spherical' 
metric on the state space. In other words, the fleld generated by the tangent vectors of the 
homogeneous coordinates Q = djtp"" is given by the gradient 

c; = -^v^if, , (4) 

where Hj is a globally deflned function on i?P", given by the expectation value Wj^^ip^ipa 
of the Hamiltonian operator Hjab in the state and V" = g°'^'Vb denotes the gradient 
operator with respect to the homogeneous coordinates. This implies that if we consider the 
coordinate transformation induced by changing each of the coupling constants 6^ according 
to the differential equation (^, then in local coordinates the vector fleld = djip"' for each 
given j, in local coordinates, satisfles the general equation for projective transformations: 

+ RlaC = ^IJb) , (5) 

where 0a is a covector given by (pa = {C\);a, Rlad is the curvature tensor, and the expression 
Xq denotes the covariant derivative of X with respect to the natural spherical metric. 

Having in mind the fairly general results noted above on the symmetry associated with 
the parameter development of the equilibrium states ip"", we would now like to specialise 
further and consider the flow induced by real-space renormalisation transformations. In 
particular, we now view the coupling constants 9^{t) themselves as dependent on a single 
scale parameter t. As a consequence, we may view the gradient flow in (Q) as a one-parameter 
family of flows, corresponding to each given scale. In this case, it is not difficult to show 
that the induced flow due to a change of scale is not quite a Hamiltonian gradient flow in 
an ordinary sense, and we require some appropriate modiflcations. 

To develop this further, we simply apply the chain rule in the equation (H) for the 
parameter development of the state vector. Hence, if we write ip"" for dip°'/dt, we have 
tp"" = —^9^Hjfjip^. Furthermore, by noticing that V^Hj = 2Hjjjip'^, we can write the deflning 
equation for the integral curve induced by renormalisation transformations in the form 

^^p'' 1 1 

^ = --V^K + -H^V^pi^ (6) 

where the potential function $ is given by the contraction $ = P^Hi of the Hamiltonian 
function with the RG /3-functions, given by = 39^ /dt. One observes that this integral 
curve does not in general correspond to a gradient flow of the form (^). However, we note 
that the potential function $ is scale dependent, i.e. $ = and is a global function on 
i^P" xR_|_. If we deflne the restriction of this function to the state space at a constant but 
arbitrary scale t = to as $4^, then the generating vector fleld for the flow, restricted to the 
given constant scale, takes the form 



dt 



= -^V«$,„ . (7) 

to 
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From this expression, we can make the following observation, namely, the flow induced by 
the renormalisation transformations on the state space manifold RP"^ induces, at each given 



scale, a projective vector fleld [^. However, the integral curves associated with such a 
flow cannot be obtained from any global projective vector fleld, since the flow changes with 
respect to the given scale characterised by the variable t. 

The situation just described is quite analogous to the mechanics of a time dependent 
Hamiltonian |T3]. In this case, the Hamiltonian 'gradient' flow gives rise to, at each instant 



of time, a Killing vector fleld. However, the integral curve along any given time evolution 
does not correspond to a global isometry of the manifold, since the Killing fleld changes in 
time. 

The general property of the induced flow in the state space that we have observed is 
worth bearing in mind when we specialise to study the flow restricted to the parameter 
space itself. In particular, since the parameter space manifold is a subspace of the state- 
space, it is natural to consider how the above structure might naturally project down onto 
the parameter space. In order to address this question we now study an explicit example, 
namely, the one-dimensional Ising model in an applied magnetic fleld. In this case the space 
of coupling constants is two-dimensional, a submanifold of the possibly inflnite dimensional 
state space, endowed with the natural Fisher-Rao metric. 



III. PARAMETER SPACE GEOMETRY FOR THE ISING CHAIN 

Consider the dimensionless Hamiltonian for the one-dimensional Ising model, given by 

N N 

H = -ir^a,a,+i-/i^(T, . (8) 

i=l 1=1 

The probability distribution for the energy H is given, as usual, by the Gibbs measure 
p{H) = exp{—H)/Z, where Z = Tr[exp(— is the partition function. Here, the symbol 
Tr denotes summation over the energy levels corresponding to all possible conflgurations. 

The totality of the space of such probability distributions {p{H)} can be mapped into a 
real Hilbert space H, as discussed earlier, by taking the square root, whereupon we obtain a 
state vector for the system iIj°'{H), which satisfles the normalisation condition gab'ip""4^'^ = 1, 
determined by the Hilbert space metric. This implies that ip"" is an element on the unit 
sphere S in Ti. The relevant parameter space submanifold Ai on the sphere is determined 
by iIj"'{6) where {6^} = [K, h) are local parameters. 



A Riemannian metric on A4, induced by the spherical geometry of S ||T0|], is given by 
the Fisher-Rao metric which, in local coordinates, is expressed as 

G,, = 4gatd,rdjij\ (9) 

where di = d/d9^ denotes differentiation with respect to the coordinates {9^,9"^) = {K,h). 
In particular, when the distribution is of Gibbs type, the Fisher-Rao metric reduces to the 
form Gij = —didj In Z, which does not explicitly involve the trace operation. 

In the case of an A^-spin Ising chain, the components of the Fisher-Rao metric can be 
calculated explicitly and are given by 
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(10) 



where r] = ysinh /i + e""^-^, and we use the notation c = cosh/i, s = sinh/i, and 
C2 = cosh2/i. The components correspond to the result obtained when we take the 
thermodynamic hmit N —>■ oo directly for the ratio Gij = limN^oodidj In Z/N. The result- 
ing expression for the metric, obtained by taking this limit, is given by M 



Gn 


= {8s 


Gl2 


= 2se 


G22 


= ce~ 



rf) + 4e 



'/ ) 



c + 77)" 



(11) 



As we have noted above, consideration of the Riemannian structure in the thermodynamic 
limit requires care, as the standard definition of the Fisher-Rao metric then requires a 
suitable integration measure over an infinite dimensional space. That is to say, the Gibbs 
measure may not be well defined as such when we have an infinite number of degrees of 
freedom. Although it appears that one may take an appropriate limit in the present case, 
some caution is required in interpreting the geometry. Indeed, as we shall discuss shortly, 
the geometry for finite N differs markedly from that observed in the thermodynamic limit. 

We first note that the components of the Christoffel symbols arising from the metric in 
this case are given by the simple expression F*^ = ^G^'-diGjk, and hence the Riemann tensor 
reduces to the form 



-pi -pm 

ml jk 



-pi -pm 

mk jl 



(12) 



In particular, in two dimensions it suffices to consider just the scalar curvature, given by 
R = 2i?i2i2/det(Gjj). Hence, by evaluating the nonzero component of the Riemann tensor 
in the thermodynamic limit, we find -R1212 = 2[e^^ {cosh h + r])]~^, from which it follows 
R that 



R 



1 + 7] ^ cosh h 



(13) 



which is always positive. 

However, if we consider the parameter space geometry for finite N using the expressions 
(p!0|), then we observe that when the size of the system N is below certain critical values 
Nc{K,h), the scalar curvature is negative. This is illustrated in Fig. 1 where the scalar 
curvature is plotted against the system size for the parameter values K = h = 1. As 
A^ increases the curvature tends asymptotically to the thermodynamic value R, while for 
N < Nc {Nc = 7 in this case) the curvature is negative, indicating a radically different 
parameter space geometry. 
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FIG. 1. A plot of the Ricci scalar R{N) at the pomt K = h = 1. The value in the thermody- 
namic limit R ~ 2.30 is approached as becomes large. 

Nevertheless, since we are only interested in the limit of large N, and as this limit 
appears to be well defined, we shall generally work henceforth with the Riemannian structure 
obtained after taking thermodynamic limit. 

We conclude this section with a brief discussion of the physical significance of the Ricci 
curvature. This is a natural question to address since the curvature can be expressed p| 
in terms of a combination of higher order central moments of the Hamiltonian H, while 
geometrically it is an invariant quantity with respect to reparametrisation. The expression 
based on the central moments of H allows a study of the scaling behaviour of R in the vicinity 
of the critical points for those models exhibiting phase transitions, from which Ruppeiner 
P] conjectured that the curvature is a measure of the correlation volume of the system (see 



In the present case we can compare R directly with the correlation length ^ of the Ising 
chain. If we write cot(20) = e"^^ smh.{h), then the correlation function C{r) = (cTjaj+r) — 
{o'i){o'i+r) can be expressed, after taking the limit — * cxd, in the simple form 



where X± = e^{c±ri) are the eigenvalues of the transfer matrix V discussed below in Section 
4. Therefore, for the correlation length ^ we find that 



where we have used (^) to determine the exact relation between ^ and the Ricci scalar R. 
Using standard techniques for inverting a power series, this relation may be expressed in the 
form. 



also i). 




(14) 




(15) 



k=0 



oo 




R 3i?2 



1 1 




(16) 
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where the coefficients are given by = — I]p=i ^ri^fe-p with the initial value Cq = R/2. 
From this relation we observe that near the critical point ^ ^ R/2 and thus R/2 provides 
a good quantitative measure of the correlation length. Such a connection breaks down, 
however, in any regime where i? ~ 1. 

Having discussed the Riemannian geometry associated with the parameter space of the 
ID Ising model, we shall, in the next section, determine the exact diffeomorphism of this 
manifold induced by an infinitesimal rescaling of the lattice spacing. Subsequently, we shall 
then consider the symmetry properties of this mapping. 



IV. CALCULATION OF THE RG /3 FUNCTIONS 

The standard approach to real-space renormalisation is implemented by means of a 
discrete scaling transformation, known as block-spin decimation. In the case of the ID Ising 
chain this transformation can be performed exactly. However, such transformations are not 
advantageous in seeking differentiable structures on the parameter space, such as we are 
investigating here. Therefore, in order to extract a generating vector field induced by these 
transformations, an appropriate form of analytic continuation is required. In this section we 
present a relatively simple technique based on the transfer matrix for achieving this aim. 

We consider first a discrete decimation of the A^-spin system via the blocking of lattice 
sites by an integer factor I at each stage of the decimation process. Thus, if we denote by a 
the lattice spacing between the spins, then at each iteration the lattice size increases by a 
factor of I, that is, a —>■ la. In addition, we assume a periodical boundary condition on the 
lattice. The renormalised parameters K' and h', after rescaling, are determined implicitly 
by the relation 

Zr,/i{K',h') = fZM{K,h), (17) 

where /, a function of K and h, is an overall scaling factor. Representing the partition 
function by Zn = Tt{V^) in terms of the transfer matrix V, given by 

V = y„ . (18) 



we can reexpress the RG relation (}[7\ ) in terms of individual configurations as 

V{K',h') = f/''V'{K,h) . (19) 
In particular, we see that the following relations 



" {v%2 ' - {v^)h ^ ^ 



hold. Therefore, in order to extract the recursion relations we require the transfer matrix 
raised to an arbitrary power I. This is achieved by the use of a similarity transformation S 
to diagonalise the transfer matrix V. Since the eigenvalues of the matrix V are given by 

A± = e^{coshh±ri) , (21) 
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if we denote the diagonalisation of V by A, i.e. V = SAS ^, we find 

= 5^^+ (22) 

Recalling the definition cot (20) = e^^ sinh h, this similarity transformation can be expressed 
in a simple form: 

S = (,\~'';''^] ■ (23) 



Thus, from the expressions in (|20|) , the recursion relations for the rescaled coupling constants 
can be obtained as 



g2/x' ^ (c + r])' + (c-r/)' tanV 

{c-r]y + {c + riY tan'^ ^ ' ^ ^ 



''''' ^ ^ ^ [ic+vy-{c-vy? ■ ^^^^ 

Here, for clarity, we recall the notation s = sinh h and c = cosh h. One readily verifies that 
these expressions reproduce the standard recursion relations for 2-spin blocking if we set 



/ = 2, and also the recursion relation for arbitrary I, when h = 0, as discussed in [|I4 . 

In the present context we assume that we can analytically continue these results to 
arbitrary / G R+. It is convenient therefore to proceed by setting Z = 1 + e, where e <^ 1. 
Then in expressions (^) and (^) the parameters now have a full lattice spacing dependence 
of the form 9 = 9{a/ao) and 9' = 9'{a/aQ + ea/ao) where 9 = {K,h) and ao is the short 
distance cutoff scale. We then perform a Taylor series expansion of the left hand side of 
(p^) and (p5| ) to 0(e). The right hand side can be expanded to the same order using the 
standard relation x"^ ~ 1 + e Inx + O(e^). Equating terms of 0(e) we obtain the /3-functions 
for the parameters K and h. In particular, it is convenient to write the lattice spacing as 
a/ Go = e*, so that the expressions = ad9^/da take the form, 

_ e~^^ sinh 2K cosh h ^ cosh h — rj 

dt 2r] cosh h + r] 

— e"^-'^ sinh 2K sinh h ^ cosh h + rj 

dt T] cosh h — 7] 

Note that and (3^ have a simple proportionality relation = —2tanh.hf3^. From these 
expressions one may readily verify various limiting cases. In particular, for arbitrary h there 
is a line of trivial fixed points at A' = 0. In the limit h —>■ 0, the first component reduces 
to the familiar form /5^q = | sinh2i^ ln(tanhi^'), while if we work to 0{h) we also find that 
(jh = -2/i/3f=o + 0{h^), as discussed in [g. 
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FIG. 2. A field plot of the differential RG generating vector field ^* in the {K, h) plane, asso- 
ciated with the one-dimensional Ising model. 

With these generating functions at hand we can view the infinitesimal RG transforma- 
tions as inducing the following diffeomorphism on the parameter space manifold, 



6' 



6' + eC 



(28) 



where e is an infinitesimal parameter and the generating vector field vec^ is given by the 
RG /3-functions, 



d_ 



(29) 



The behaviour of this vector field is displayed in Fig. 2, indicating the flow of the RG tra- 
jectories away from the critical point at K = oo towards the line of stable high temperature 
fixed points at K = 0. 

In the next section we shall study the geometrical properties of this vector field in various 
limiting cases. Before proceeding to do so, we note that certain physical features of the model 
contained in the vector field ^* are only made explicit when the underlying Riemannian 
structure is taken into account. To this end, we consider the conventional definition of the 
fixed points, as identified with the zeros of the /5-functions. This definition is not entirely 
satisfactory, since there are cases where the /3-function does not vanish at the critical point, 
such as the example considered here. While the /3-functions vanish at one of the fixed points, 
K = 0, they tend to constant values for any given value of h when K ^ oo. 

Clearly, if the underlying geometry were Euclidean, then this implies that {K = oo,h = 
0) cannot be regarded as a proper fixed point. However, the point here is that the parameter 
space is a curved manifold endowed with the Fisher-Rao metric. Therefore, any fixed point 
of a given vector field generating a diffeomorphism on the manifold, in general, can be (cf. 



|15|) identified with the zeros of the 'velocity function' defined by 



(30) 



The behaviour of the velocity function in the case of the ID Ising model is illustrated in 
Fig. 3. It is clear that the proper velocity of the fiow vanishes at both fixed points. An 
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alternative interpretation of this point is available in two dimensions where the existence of 
the Zamolodchikov C-function allows the (3 function to be constructed as /?* = G^WjC near 
the fixed points |T^. When the renormalization group generates a gradient flow of this form 
it is then clear that the fundamental covector j3j 
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djC also encodes the correct fixed 



point structure of the flow. 




FIG. 3. A plot of the velocity function v{K, h) over the parameter space of the one dimensional 
Ising model. The velocity of the flow vanishes for an arbitrary value of h at the high temperature 
fixed points at IT = 0, and also at the critical point given by ii' = oo and /i = 0. 



V. THE SYMMETRY OF RG DIFFEOMORPHISMS 

As we have discussed in Section 2, the RG generating vector field in the projective 
state space may be interpreted as generating a projective automorphism of the fixed scale 
state-space manifold. On the other hand, if we consider the global flow induced by the 
RG, then there seems to be no generic symmetry associated with the flow that can be 
expressed in terms of conventional diffeomorphisms in Riemannian geometry. This fact then 
leads us to conjecture that such a symmetry should also be absent for RG flow viewed in 
its normal setting of the parameter space submanifold. In this case the restriction of the 



natural state-space metric to the parameter space can be shown |jTT| to coincide with the 
Fisher-Rao metric which, as we have discussed in Section 3, can be explicitly evaluated 
in local coordinates which are parameters of the Hamiltonian. The components of the 
generating vector field for the full scale-dependent flow then reduce to the /3-functions for 
these parameters. Thus we are led to enquire whether the vector field ^* is a projective vector 
field on the parameter space manifold. Indeed, this is also a natural question on general 
grounds from the standard theory of Riemannian geometry, since projective automorphisms 
are the most general diffeomorphisms preserving geodesies on a manifold, and include affine 
and Killing transformations as special cases. 

In fact, as we shall observe below, by explicit analysis of a counter example — the one 
dimensional Ising model — the fiow induced by the RG in general does not possess any 
standard global symmetry. In particular, the exact generating vector field for RG fiow in 
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the ID Ising model, is not a projective vector field. This is of course consistent with the 
observations in state space where there is an explicit correction term. 

This statement is verified in the case of the ID Ising model in a straightforward manner 
since for ^\ evaluated in Section 4, to be a projective vector field it must satisfy the relation 

(I), 

4, + 4,e^ = 5(';A), (31) 

for some covariant vector (pi. In two dimensions, the components of this equation are ex- 
plicitly given by 



^;2;2 + Rl221^^ = ^fi-i + G^^-R2112^^ = 0, 

+ G'^R2U2e = 201 ef2;2 + G^'Ri22ie = 202 

^;\;2 + ^^^-^2121^^ = 02 ^;2;1 + G'^^ Ri2l2^^ 

^;2;1 + ^^^-^1212^^ = 02 ^fl;2 + G'^^-R2121^^ 



(32) 

'1, 
>i, 



and by direct substitution of the expressions for the metric, curvature, and the vector field 
C,^ , we find that the equation is not satisfied globally on the parameter space for any choice 
of 0j for the geometry at finite N and also in the thermodynamic limit. Thus is not a 
projective vector field on the parameter space manifold. Since we have knowledge of the 
exact geometric quantities, and the RG /3-functions, in this case the ID Ising model in an 
external field serves as a counter-example to any conjecture of the general validity of such a 
symmetry, as stated above. 

It is interesting to note that as we are dealing in this case with a two-dimensional pa- 
rameter space, this allows additional control through the fact that one may coordinatise the 
surface with a complex Riemann coordinate and consider general conformal transformations. 
However, as we wish to consider this model purely as a tractable example of higher dimen- 
sional (generically infinite dimensional) cases we shall avoid for the moment any aspects 
which appear particular to two dimensions. 

Nevertheless, we can consider various limiting cases by restricting attention to local 
regions of the parameter space. The most relevant region is of course the vicinity of the 
critical point. Therefore, it would be convenient to linearise the vector field around K = oo 
and h = 0. However, as we notice from the relations p6| ) and (pTf), the generic dependence 
on the temperature parameter T = is of the form exp(±2/T), implying that T = is an 
essential singularity. Thus, following standard practice for the study of critical exponents, 
it is convenient to introduce a new variable r = exp{—2K), and consider a linearisation of 
the RG transformation in the variables {T,h) about (0,0). Assuming that r and h are of 
0(e), and expanding equations (^6]) and (^) up to 0(e''), we obtain 

3 = h + \ \ + 0(en , (33) 

^ 2 V 3 6 ; V 15 15 9oy ^ ^ ' ^ ^ 

= /.-(^)+0(6^). (34) 

Note that various geometric quantities, e.g., the curvature, remain singular at the critical 
point in these variables. However, in terms of r and h we can construct a Laurent expansion 
and consider the structure of the lowest order terms. 



13 



Rather than considering general projective automorphisms in the neighbourhood of the 
critical point, it is now convenient to analyse particular transformations which have had some 
attention in the literature. In particular, we consider the degree to which RG flow violates the 
geodesic equation, and whether there are any special flows which are exact geodesies of the 
metric. We also investigate the extent to which ^* may be regarded as a conformal Killing 
vector near the critical point. Conformal Killing transformations, preserving the angles 
between vectors, are not a subset of projective automorphisms but one readily verifies that, 
for the ID Ising model, the differential RG does not generate such a symmetry globally. 

We first analyse the relationship between RG diffeomorphisms and geodesic flow on 
the parameter manifold. In fact, using the standard discrete decimation procedure, Dolan 
suggested that the induced flow might be a geodesic along the line h = 0. Although 
an explicit solution for the geodesic equations has not been obtained, we can nevertheless 
determine in which regimes the vector field ^* induced by the infinitesimal scale change 
satisfies the geodesic equation ^-^ Vj^' — A^'' = 0, where Vj denotes the covariant derivative 
compatible with the underlying Fisher- Rao metric Gij. The function A in this expression 
vanishes only if the parameter t = Ina/ao is an affine parameter. 

Numerical analysis indicates that along the lines h = and K = the flow is an 
exact geodesic flow, and that the deviation becomes linear in h, as K increases. Hence our 
results confirm the expectation in [Q. In order to see this explicitly, we expand the geodesic 
equations in r and h to obtain 

eV,e^ -Ae = + + ^ + 0{e') , (35) 

e V,e' -A^^ = + (1 - A)h + 0{e^) . (36) 

The common factor of A is given hj A = 1/2 + 0(e), and the deviation from geodesic flow 
is clearly observed to be linear in h with a positive coefficient of 1/2 for each component. 

We turn now to an investigation of whether corresponds to a conformal Killing vector 
field in particular regions of the parameter space. Justification for the conjecture that ^* may 
have such a structure at least locally on the manifold in the neighbourhood of the critical 
point is provided by the analysis of Diosi et al. who have shown that one recovers the 
standard critical exponents of a system such as the ID Ising model under the assumption 
that linearised around the critical point is a conformal Killing vector field. In other words, 
the components of the metric tensor Gij should satisfy C^G — dG = where denotes the 
Lie derivative with respect to the vector field ^, and d is the spacetime dimension. 

As we have noted earlier, explicit evaluation of this equation in the present example 
indicates that the symmetry does not hold globally on the parameter space. However, is 
indeed a conformal Killing field to a good approximation near the critical point K = oo, 
h = 0. Indeed significant deviations are only seen as one moves close to the high temperature 
fixed points where K = 0. 

The approximate symmetry near the critical point may be shown more quantitatively 
by expanding the conformal Killing equation in a Laurent expansion in the variables r and 
h. For example, the (1, 1) and (2, 2) components take the form 

{C,G)n - ^Gn = + f - 8(2 - df-^^^ + 0{h')] + O(e^) , (37) 
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{C,Gh, - dG,, = + - l^^^m^ + Oih')) + 0{e^) , (38) 

and thus ^* in fact corresponds to a Killing vector field up to 0(r), and to a conformal 
Killing vector field with d = 1 up to O(r^) for the (1, 1) component and up to 0{r'^/h) for 
the (2, 2) component. We note that this expansion is not a strict linearisation of the RG 
equations about the critical point, due to their singular structure. Nevertheless these results 
verify, in this example, the claim of Diosi et al. explicitly. 

VI. DISCUSSION 

We have presented an analysis of the Riemannian parameter space geometry of the ID 
Ising model, and the characteristics of renormalisation group trajectories on this manifold. 
While the generator of this transformation is, up to higher order corrections, a conformal 
Killing vector field near the critical point, there is apparently no such interpretation globally 
on the manifold. This is consistent with the general structure observed for RG flow in the 
state space of the system. 

It should be pointed out, however, that in other models of interest such symmetries 
may exist due to particular properties of the theory itself or the individual parameters. 
The question of the existence of other independent parameter space symmetries is also of 
importance as these may provide additional, generically nonperturbative, constraints on the 
/3-functions of the theory. Discrete parameter space symmetries are of this kind, while 
in the present case we observe that there is a class of diffeomorphisms of the parameter 
space manifold which preserve the RG flow in the sense that, for the generator X of such a 
transformation, CxV = where r] = rjidO^ is the 1-form dual to the vector field One may 
verify that 7] has the form rj = (?7i,0) and thus any diffeomorphism generated by a vector 
field of the form X = (0, X2) will preserve the RG flow in the sense described above. Such a 
structure is clearly quite model dependent. Nevertheless, the presence of a commuting flow 
may well have some more general validity. 
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